Magnonlike Dispersion of Spin Resonance in Ni-doped BaFe 2 As 2 



CO 

i— I 

o 

CN 

u 

< 

CN 



o 

o 
i 



c 
o 
o 



> 

in 

O 



s 



M. G. Kim, 1 '* G. S. Tucker, 1 D. K. Pratt, 1 S. Ran, 1 A. Thaler, 1 A. D. Christianson, 2 K. Marty, 2 S. Calder, 2 

A. Podlesnyak, 2 S. L. Bud'ko, 1 P. C. Canfield, 1 A. Kreyssig, 1 A. I. Goldman, 1 and R. J. McQueeney 1 ' 1 

'Ames Laboratory and Department of Physics and Astronomy, Iowa State University, Ames, IA, 50011, USA 
2 Oak Ridge National Laboratory, Oak Ridge, TN, 37831, 
(Dated: April 25, 2013) 



USA 



Inelastic neutron scattering measurements on Ba(Feo.963Nio.o37)2As2 manifest a neutron spin 
resonance in the superconducting state with anisotropic dispersion within the Fe layer. Whereas 
the resonance is sharply peaked at Q APM along the orthorhombic a axis, the resonance disperses 
upwards away from Qafm a l° n g the b axis. In contrast to the downward dispersing resonance and 
hour-glass shape of the spin excitations in superconducting cuprates, the resonance in electron-doped 
BaFe2As2 compounds possesses a magnon-like upwards dispersion. 

PACS numbers: 74.70.Xa, 74.20.Mn, 78.70.Nx 



A spin resonance, observed in inelastic neutron scatter- 
ing measurements, appears in superconducting materials 
that do not possess a simple s-wave gap symmetry. [1, 2] 
Consequently, the spin resonance is considered to be a 
hallmark of unconventional superconductivity and high- 
lights the important relationship between antiferromag- 
netic spin fluctuations and superconductivity. The spin 
resonance has an intimate connection to the super- 
conducting state. [1-5] The inelastic neutron scattering 
(INS) signature of the resonance is a gapping of normal 
state spin fluctuations at the antiferromagnetic (AFM) 
wavevector, Qafm; f° r energies well below the super- 
conducting gap 2 A, and an enhancement of magnetic 
spectral weight at an energy fio < 2A. The energy 
of the resonance at Qafm i s found to be proportional 
to the superconducting (SC) transition temperature T c 
{no/k B T c « 4 - 6) or 2A (ft /2A « 0.6) for a variety 
of unconventional superconductors. [6, 7] The dispersion 
of the spin resonance, f2 q (where q = Q — Qafm), is a 
measure of the resonance energy away from Qafm ■ 

In the case of the single-band cuprates, the disper- 
sion of spin fluctuations forms a characteristic hour-glass 
shape below T c where the resonance appears at the neck 
of the hour-glass with energy Hq, as shown in Fig. la. [3, 4] 
The downward dispersion (u < fi ) is associated with 
the resonance itself, 17 q , whereas the upward dispersion 
(w > fio) arises from normal state spin waves. [8] In the 
cuprates, the downward resonance dispersion is deter- 
mined by nodes in the d-wave SC gap at q = q no dc where 




A qnod< 



= 0. 



Although the spin resonance has been clearly observed 
at Qafm m multi-band iron arsenide compounds stud- 
ied to date, no dispersion of the resonance within the 
Fe layers (in the ab plane) has yet been reported and, 
complicating matters, band structure calculations pre- 
dict both an upward [9] and downward [10] dispersion 
depending on the details of the bands and the symmetry 
of the superconducting order parameter. In this Letter, 
we use inelastic neutron scattering to show that the res- 
onance in Ba(Fe .963Ni .o37)2As2 (T c = 17 K) disperses 
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Figure 1: (a) Spin fluctuations in the cuprates showing the 
hour-glass shape of the spin excitations. The downward dis- 
persing part (orange line) is the resonance mode for a d-wave 
gap and the upward part (black line) is the competing AFM 
spin fluctuations, (b) Spin fluctuations in our Ni-doped iron 
arsenide superconductor where both the resonance and com- 
peting AFM fluctuations disperse upwards. 



upwards. Using a simple mean-field approach with the 
assumption of s ± order parameter, we show that the de- 
tails of the dispersion, such as the resonance velocity, are 
determined by the normal state spin fluctuations. 

The sample studied is a single crystal of 
Ba(Fei_ x Ni x ) 2 As2 with x — 0.037 weighing 436 mg. 
Rather than using co-aligned single crystals with a 
larger total mass, we chose to measure a single spec- 
imen of high crystallinity (mosaic width < 0.44°) to 
minimize effects of sample mosaic on the q-dependence 
of the spin fluctuations. Our preliminary neutron and 
x-ray scattering measurements show that the sample 
orders into an incommensurate AFM structure [11] 
below Tn = 26 K with a tetragonal-orthorhombic 
transition below Tg — 29 K. The superconducting 
transition temperature is T c = 17 K. Additional details 
of crystal growth and characterization are described 
elsewhere. [12] The sample was mounted in a closed- 



cycle refrigerator for temperature dependent studies 
on the HB3 neutron spectrometer at the High Flux 
Isotope Reactor at Oak Ridge National Laboratory 
using pyrolitic graphite (PG) (0, 0, 2) reflections for 
both monochromator and analyzer and a configuration 
with E f = 14.7 meV, 48' - 60' - 80' - 120' collimation 
and two PG filters before the analyzer. We define 
Q = (H,K,L) = ^-Hi+ ?Z-Kj + ^-Lk where the 
orthorhombic lattice constants, a ~ b = 5.6 A and 
c = 13 A were determined at T = 15 K. The crystal 
was aligned in the (1,0,1) — (0,1,0) scattering plane. 
Raw data were converted into the imaginary part of 
the magnetic susceptibility, x"(Q. w ) after subtracting 
estimates of the non-magnetic background and correct- 
ing for the Bose population factor. In addition, we 
performed inelastic neutron scattering measurements on 
the Cold Neutron Chopper Spectrometer (CNCS) at 
the Spallation Neutron Source at Oak Ridge National 
Laboratory using an incident energy of 10 meV. Data 
were accumulated in the [H, K, 1] plane by rotating the 
crystal around a vertical (0, 1, 0)-axis and performing a 
series of exposures. 

Figure 2a shows the energy dependence of the imag- 
inary part of the magnetic susceptibility, x"(Q AFM ,o;), 
for Ba(Feo.963Ni .o37)2As2 in the normal state at T = 
20 K and the superconducting state at T — 4 K, where 
Qafm = (1) 0, 1) in orthorhombic coordinates. The neu- 
tron spin resonance peaks at an energy of fl Q s=s 6 meV, 
as made clearer by the difference plot in Fig. 2b showing 
the resonance enhancement. At Qafm. the resonance ex- 
tends up to approximately 10 meV, above which the sus- 
ceptibilities of the normal and SC states are equivalent. 
Measurements of the energy spectrum at positions offset 
from Qafm m the transverse direction [Q = (1, —0.1, 1) 
and (1, —0.15, 1)] show that the center of the resonance 
intensity shifts up to higher energy (Figs. 2c— f). In ad- 
dition, we observe that the resonance spectral weight ex- 
tends to at least 14 meV. Thus, despite the very broad 
lineshapes, measurements of the susceptibility away from 
Qafm show that the spectral weight of the resonance has 
moved to higher energy, i.e. the resonance is dispersing 
upwards, unlike the cuprates. 

The observation of resonance dispersion is confirmed 
by a series of constant-energy Q-scans in the normal and 
SC states in the vicinity of Qafm = (1j 0' 1)> focusing on 
the transverse direction Q = (1,K,1) [q = (0,K, 0)] as 
shown in Figs. 3a and b. Here, the resonance is clearly 
identified by the enhanced susceptibility at Qafm an d 
to s» Slo = 6 meV. As the energy is increased, the spectral 
weight progressively moves away from Qafm an d weak- 
ens. Plots of the difference between the susceptibility 
in the normal and superconducting states more clearly 
show that the resonance, fi q , disperses upwards away 
from Qafm (Fig- 3b). In contrast, along the longitu- 
dinal (H, 0, H ) direction in our scattering geometry, we 
find that only a weak resonance enhancement remains at 




3 6 9 12 1 
co(meV) 



3 6 9 12 1! 
(o (meV) 



Figure 2: The imaginary part of the magnetic suscepti- 
bility measured on HB3 at (a) Qafm — (1, 0, 1), (c) Q = 
(1, —0.1, 1), and (e) Q = (1, —0.15, 1) in the normal state at 
T = 20 K (filled circles) and the superconducting state at 
T — 4 K (open circles). Red lines are model calculations of 
the normal state susceptibility assuming the NAFL spin fluc- 
tuations, as described in the text. The difference between 
superconducting and normal state susceptibilities is shown at 
(b) Qafm. (d) Q = (1,-0-1,1), and (f) Q = (1,-0.15,1). 
Black lines are guide to the eyes and arrows indicate the res- 
onance peak obtained from Gaussian fit. 



Qafm a ^ 9 meV and no resonance enhancement is ob- 
served at 12 meV, i.e. there is surprisingly no indication 
of dispersion in the longitudinal direction (Figs. 3e and 
f). This observation is consistent with early measure- 
ments of the resonance in both Ni (ref. 13) and Co-doped 
(ref. 14-16) BaFe 2 As2 that find a sharply defined reso- 
nance in the H (orthorhombic a) direction. Thus, while 
the resonance dispersion is observed in the transverse di- 
rection, it is apparently too steep to observe a splitting in 
the longitudinal direction due to the limited instrumen- 
tal resolution (we will return to this point below). As 
discussed below, the Q-space anisotropy of the spin reso- 
nance in the SC state may be associated with the normal 
state spin fluctuations within the Fe layers that possess 
a two-fold (elliptical) anisotropy (Fig. 3e). [17-20] 

We now turn to a discussion of the spin dynamics in 
the normal and SC states, which leads to a straight- 
forward interpretation of the resonance dispersion and 
its anisotropy. Despite the presence of weak AFM or- 
dering, the normal state spin dynamics of the under- 
doped BaFe2As2 compounds at low energies can be un- 
derstood using the model of nearly AFM Fermi liquids 
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Figure 3: (a) The imaginary part of the magnetic susceptibility measured on HB3 at various energy transfers along the [0, K, In- 
direction transverse to Qafm = (L u i 1) i n the normal state at T = 20 K (filled circles) and the superconducting state at T — 4 K 
(open circles), (b) The difference of superconducting minus normal state susceptibilities from (a), (c) The susceptibilities along 
the [H, 0, H] -direction longitudinal to Q AFM in the normal (filled circles) and superconducting states (open circles) and (d) 
the difference. Red and blue lines are model calculations of the normal state susceptibility and superconducting resonance, 
respectively, as described in the text, (e) CNCS measurements of the normal state spin fluctuations in the (H, K, 1) plane at 
w = 6 meV and T = 25 K. (f) Fit of the CNCS data to a model of NAFL, as described in the text, (g) Contour plot showing the 
difference of the magnetic susceptibilities in the superconducting and normal states measured on HB3 depicting the dispersion 
of the resonance in the transverse direction. Symbols correspond to fitted values of the resonance maxima and the line is a fit 
to the magnonlike dispersion given in Eq. 4. The vertical error bars indicate the energy resolution of HB3 and horizontal error 
bars indicate the full-width-half-maximum of fits. 



(NAFL). [17, 21] The imaginary part of the dynamical 
susceptibility, as measured by INS, can be written as in 
ref. 21, 



x£(q>w) 



xo^r q 



r q [i + q g 2 



(i) 



where £ q is the AFM correlation length, T q is the relax- 
ation width due to the decay of spin waves into electron- 
hole pairs (Landau damping) , and the subscript q allows 
for anisotropy in these quantities. 

The anisotropy of the normal state susceptibility 
has been studied in detail. [19] The two-fold in-plane 
anisotropy of the magnetic correlation length is most im- 
portant for the subsequent discussion, 



Li = £ W + V{Qx - Qy)\ 



(2) 



and we assume in Eq. 1 that the damping T q f=a Tq 
is isotropic and constant throughout the Brillouin zone. 



Since the resonance dispersion observed along the L (c) 
direction in AFe2As2 {A — Ca, Sr, and Ba) based su- 
perconductors is weak, we assume q z gives a negligible 
effect. Typical fits to the normal state data in terms 
of the NAFL convolved with the instrumental resolution 
are shown in Figs. 2, 3a, 3c, and 3f with Tq fixed to 10 
meV. [19] We obtain the following parameters; £ = 9.5(4) 
A and r\ — 0.5(1) corresponding to different correlation 
lengths 11.8(7) A and 6.5(3) A in the H (orthorhombic 
a) and K (b) directions within the Fe layers, respectively. 

The neutron spin resonance is reasonably well un- 
derstood as an excitonic bound state in an itinerant 
antifcrromagnet.[l, 2] In this picture, fermions form sin- 
glet Cooper pairs and electron-hole (singlet-triplet) exci- 
tations appear with a threshold energy | Ak|+| Ak+Q AFM |. 
These electron-hole excitations will form a bound-state 
inside the SC gap (an exciton) as a consequence of inter- 
actions already present in the normal state. The spec- 



tral features of the bound-state can be understood within 
a mean-field or random-phase approximation (RPA) ap- 
proach to the interactions, as described below. 

Following closely the arguments from Ref. 22, the RPA 
analysis of the susceptibility in the superconducting state 
leads to an equation for the dispersion of the resonance 
mode, 



tii 



A q r q (i 



2„2\ 



Qz 



(3) 



where 2A q = |A k | + |A k+ q AFM+q | is the fermion gap at 
two points on the Fermi surface connected by q + Qafm 
and fio = V^oTo is the resonance energy at Qafm 
(q = 0). Note that for a single-band model with a d- 
wave gap possessing nodes at q no dc, Ao is a maximum 
and A qnodo — > at points on the Fermi surface, thereby 
forcing ^q nodo — > and resulting in the downward disper- 
sion shown in Fig. la. The simple RPA model therefore 
captures the essential features resulting in the downward 
dispersion in the cupratcs. 

Although the q-dependence of the resonance in Eq. 3 
has several contributions (such as q-dependent gap or 
damping effects), here we choose to apply the RPA ap- 
proach for the iron arsenides under the assumption that 
the q-dependence arises only from anisotropy in the nor- 
mal state spin fluctuations. Assuming the proposed s 
gap symmetry for our compound, A q = Ao for all q and 
the equation above takes the form of a gapped magnon 
with upward dispersion 



Oj = 



Q2 



"t* ^res.q^ 



(4) 



where c 2 res ^ 



ilp^q is the velocity of the resonance mode 
which becomes anisotropic itself due to the anisotropy 
in the normal state spin-spin correlation length. More 
detailed four-band RPA calculations (ref. 23) and two- 
orbital calculations of the susceptibility within the self- 
consistent fluctuation exchange (FLEX) approximation 
(ref. 9) arrive at similar conclusions regarding the up- 
ward resonance dispersion, at least at low doping con- 
centrations. With increased electron doping, other multi- 
band RPA calculations have predicted the development 
of an incommensurate normal state response. [10, 20, 23] 
In the superconducting state, these calculations predict 
that the resonance itself is also incommensurate and first 
disperses downwards to the incommensurate wavevector, 
and then upwards. Our observation highlighting the up- 
ward dispersing and commensurate response is, there- 
fore, consistent with the multi-band RPA approach for 
underdoped compositions. We note that an incommen- 
surate resonance was reported for Ba(Fe . 925^0. o7s)2As 2 
by Luo et al.[2A] at a higher composition than reported 
here. This evolution towards an incommensurate re- 
sponse would be consistent with the general trend out- 
lined above. However, the incommensurate splitting was 
only observed in the superconducting state (not in the 



normal state) and found to disperse upwards away from 
Qafm rather than downwards. 

We predict that the anisotropic resonance velocities, 
yjtll£?{\ ± 77), in the Fe layers to be 85(5) meV A and 
50(5) meV A in the H and K directions, respectively, us- 
ing values (£ = 9.5(4) A and 77 = 0.5(1)) obtained from 
the normal state spin fluctuations. These predictions 
can be compared to fits made to the total susceptibility 
(x" = Xn + Xres) after convolution with the instrumental 
resolution function, where Xres l& the resonance enhance- 
ment that obeys the dispersion relation in Eq. 4 and \n 
is given by Eq. 1. Typical fits are shown as the blue line 
in Figs. 3a and b. A contour plot of Xres obtained from 
INS data along with fitted values of the peak positions is 
shown in Fig. 3g and we arrive at 63(2) meV A for the 
observed resonance velocity in the K direction, in good 
agreement with the simple analysis above. 

The resonance velocity is to be contrasted with the 
velocity of the normal state AFM spin waves which are 
much steeper c rj 450 meV A ^> c res . Within the disor- 
dered AFM state model, dispersive features of the normal 
state AFM spin waves will appear only at much higher 
energies, u> > c£ _1 w 60 meV, as shown in Fig. lb. 
This disparity in the magnitude of the normal state and 
resonance velocities eliminates the possibility that reso- 
nance may be interpreted as damped AFM spin waves 
that sharpen up in the SC state (see Supplement for 
more information). We also note that limitations due to 
experimental resolution may prevent observation of the 
splitting in the longitudinal direction with the predicted 
velocity of 85(5) meV A as shown in Fig. 3c and d. 

In conclusion, we find that the dispersion of the neu- 
tron spin resonance in Ba(Fe .963Ni .o37)2As2 can be 
interpreted based on the assumption of extended s- 
wave superconductivity and the properties of the NAFL 
spin fluctuations in the normal state, most notably the 
anisotropic spin-spin correlation length. Such obser- 
vations strongly support an excitonic picture for the 
spin resonance. Similar reasoning may also explain the 
weak resonance dispersion observed along the direction 
perpendicular to the layers (i-direction) in ^4Fe2As2 
based superconductors. For quasi-two-dimensional nor- 
mal state spin fluctuations observed in optimally doped 
samples, the resonance will have no dispersion along L 
(appear flat) due to the vanishing correlation length. 
However, underdoped samples with long-range AFM or- 
der, such as Ba(Fci_ a; CO; C )2As2 with x = 0.04 [16] 
and 0.047[15, 25], still possess weak interlayer spin cor- 
relations. Correspondingly a finite, albeit small, L- 
dispersion is observed in these samples (with a velocity 
of ~ 6 meV A). In a similar vein, doping will also af- 
fect both the interlayer and intralayer correlation length 
and consequently the resonance dispersion is expected 
to be composition dependent. This could explain recent 
observations of a large transverse resonance splitting in 
overdoped Ba(Fei_ :E Ni a; )2As2.[24] 
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Supplement: Magnon-like dispersion of neutron spin resonance in Ni-doped BaFe 2 As 2 
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The following supplemental information covers the derivation of the dynamical magnetic susceptibility of itinerant 
antiferromagnets in the normal and superconducting states and the approximations made in the derivation. 

Normal state. Despite the presence of weak AFM ordering, the normal state spin dynamics in the undcrdoped 
BaFc2As2 compounds can be understood using a model of spin waves in a disordered AFM state. In the long- 
^vq wavelength limit, the dynamical susceptibility can be written as 1 

_4- where £ q is the AFM correlation length, c q is the AFM spin wave velocity and the subscript q allows for anisotropy 

£Sj in these quantities. The wavevector q = Q — Qafm ^ s measured relative to the antiferromagnetic wavevector Qafm- 

For a spin-density wave state, II„(q, us) = «7 q w is purely imaginary and arises from the decay of spin waves into 

electron- hole pairs (Landau damping) . Neutron scattering measures the imaginary part of the susceptibility, given by 



cd 



cd 



R at \ = xoC q r q ^ 

^ Xn[ d, «>) „* + F 2 [! + q^ _ W 2^ /C 2 ]2 ^ ) 

where T^ 1 = 7 q £q is the relaxation width due to Landau damping. We have showed previously in undcrdoped 
Ba(Fe 1 _ 2 ,Co a; )2As2 that on average c w 450 meVA and £ w 7.5 A and c^ 1 sa 60 meV. 2 So for energies relevant to 
the superconducting resonance us < 2A sa 10 meV <C c£ _1 , we can ignore the dispersive spin waves and the normal 
state susceptibility falls into the diffusive limit 



■ at \ xoC q r q w 

^ xM ^ ] -^ + ri[i + e^?- (3) 

Fig. SI (a) shows a typical form of the response of an itinerant antiferromagnet in the normal state. 

Superconducting state. In the superconducting state, a similar form for the susceptibility can be used 1 

^H 2 

> Xsc (q, W )= X oK q 2 + 9 2 -^-n sc (q, W )]- 1 . (4) 

iq However, in this case the opening of the superconducting gap eliminates the imaginary contribution of the Landau 

\Q damping to the self energy. Instead, the smallest non-vanishing term is quadratic in energy, LI sc (q, us) ~ 7 q w 2 /A q . 

■ ' Here 2A q = |Ak| + |Ak+Q APM + q | is the fermion gap at two points on the Fermi surface connected by q + Qafm- The 

f» ^ susceptibility can then be written 

X sc (q,^) ~ Xo[? q 2 + q 2 u\\ + ^)Y\ (5) 

C q ZA q 

Again using the properties of the normal state spin fluctuations (the isotropic damping, To, and the AFM correlation 
length, £) 2 and the superconducting gap magnitude (Ao), we estimate that a/AoToC 2 ~ \/(5 meV)(10 meV)(8 A) 2 ss 

60 meVA. As this is much smaller than the the AFM spin wave velocity, this shows that we can again ignore the spin 
wave term and we arrive at the following approximation for the susceptibility: 

^(q^-XoCqll + ^-^r 1 . (6) 

While this form of the susceptibility contains only a real-part and is not directly comparable to the neutron scattering 
data, the resonance condition will correspond to the vanishing of the term in square brackets in the denominator, 
leading to an equation for the dispersion of the resonance mode 

fi q = A q r q (l + £ 2 <z 2 ) (7) 




2A 




q Q £l q 2A oo 

Figure SI: The response of an itinerant antiferromagnet and the resonance dispersion for an s-wave superconductor (a) 

A typical form of the response of an itinerant antiferromagnet in the normal (red shade). The superconducting response is indicated by 
the resonance dispersion (blue line). Dispersion of the normal state AFM spin waves appear above to > cf — 1 and the SC resonance occurs 
at Q.q -C eg -1 and disperses upwards, (b) A graphical explanation of the resonance dispersion for an s-wave superconductor. The solid 
line corresponds to the quadratic behavior of rise i n the superconducting state. The SC resonance appears at u) = Qo when ilgc = £ -2 
(with q 2 = 0). At q 7^ 0, the resonance will always occur at a larger energy, f2 q . 



where at Qafm (q = 0), ^o = \ZAqTq is the resonance energy. Note that for a d-wave gap, Ao is a maximum and 
A qnodo — > at the node wavevector q no( j e resulting in a downward dispersion. For the isotropic s-wave gap of the iron 
arsenides, A q = Aq for all q and the equation above takes the form of a gapped magnon 



<X 



ni 



(8) 



where c 2 es q = Aor q £ 2 is the velocity of the resonance mode. One can refer to Fig. Sl(b) for a graphical explanation of 
the resonance dispersion for an s-wave superconductor. In this figure, we see that the resonance condition corresponds 
to the vanishing of the susceptibility in Eq. 4 (or denominator in Eq. 6) where q 2 + £~ 2 — Use = 0. The resonance 
will always occur at a larger energy for non-zero values of q when compared to q = 0. 

Anisotropy. The normal state spin fluctuations in the iron arsenides possess a two-fold anisotropy about Qafm- This 
anisotropy has been studied in detail 2,3 and can be cast in the form of anisotropies in the magnetic correlation length 
and damping, 



ty = eii 2 + vM 



al)\ 



and 



T cl = r [l + a 2 (q 2 + r lr (ql-q 2 y )} 



(9) 



(10) 



where q is defined in the coordinate system of the low temperature orthorhombic structure, i]^ is the anisotropy of the 
AFM correlation length, a is a constant allowing for momentum-dependent Landau damping, and rjr is the anisotropy 
of the Landau damping. This anisotropy of the normal state spin fluctuations will lead to anisotropy in the resonance 
dispersion transverse (T, q = q y ) and longitudinal (L, q — q x ) to Qafm = (1,0, 1). The resonance velocities in these 
two directions are 



A r K 2 (l-^)+a 2 (l-r ?r ) 



(11) 



and 



c 2 reStL = A r K 2 (i + Vi) + a 2 (l + rtr)] 



(12) 



Referring to fits to the normal state spin fluctuations in a Co-substituted sample one last time 2 , we have the 
following parameters; £ =7.5 A 77^ = 0.6, Tq = 11 meV, a — 2 A and rjr = -1- We notice that, for this particular 



set of parameters, the anisotropy in the damping contributes very little to the resonance velocity and can be ignored, 
resulting in the simple forms for the velocities. 



^rea,T 

and 



A r £ 2 (l-^) (13) 



c 2 re ^ L = A Q T e(l + Vi) (14) 
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